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We consider and compare two different approaches to the fractional subdiffusion and transport 
in washboard potentials. One is based on the concept of random fractal time and is associated 
with the fractional Fokker-Planck equation. Another approach is based on the fractional general- 
ized Langevin dynamics and is associated with anti-persistent fractional Brownian motion and its 
generalizations. Profound differences between these two different approaches sharing the common 
adjective "fractional" are explained in spite of some similarities they share in the absence of a non- 
linear force. In particular, we show that the asymptotic dynamics in tilted washboard potentials 
obey two different universality classes independently of the form of potential. 

PACS numbers: 



I. INTRODUCTION 

Anomalous diffusion becomes an increasingly popular subject with the number of papers published per year growing 
fast over last twenty years after a distinct rise occurred in 1990 (cf. Fig. [T]). Since then, it spreads from such typical 
physical applications as charge transport in disordered solids and hot plasmas to biophysical applications and even 
quantitative finance^— . 

There are several sufficiently generic physical mechanisms and accompanying theoretical approaches to describe 
the complexity of anomalous transport processes. One approach is intrinsically based on the physical picture of a 
stochastic time clocki^^— It models random sojourns of a travelling particle in trapping domains of a disordered 
solid, e.g. due to energy disorder. After a random time spent in some spatially located trapping domain the particle 
jumps to a neighboring domain, or maybe farther, and such a jumping process continues in time. The jump directions 
and their lengths are not correlated from jump to jump and the next clock period is not correlated with those passed 
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FIG. 1: Number of papers on the subject "anomalous diffusion" published per year according to the ISI Web of Science, 
Thomson-Reuters. Notice the "phase transition" -like rise during 1990. 
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(semi-Markovian assumption 1 ). Such a random clock is completely characterized by the probability density of clock 
periods ip(r). In a given time interval t there will be a random number of jumps n, or stochastic clock periods 
completed. However, if the mean clock period (r) exists, the probability distribution p(n, t) of "ticking" n times 
within the observer time window t becomes for large n a very sharp function 4 around n*(t) = t/(r), as characterized 
by the relative dispersion, (5n 2 (t)) 1 / 2 / (n(t)) (see Appendix A). In the continuum medium approximation the trapping 
domains shrink to points. Then, (r) can be made arbitrarily small and n becomes quasi-continuous variable for a finite 
t. Correspondingly, the probability density of intrinsic time, r(t) = n{t)r sc , where t sc is a time-scaling parameter, 
an intrinsic clock time unit equal to the duration of time period for the regular clock, assumes a delta-function, 
p(r) = 8{t — t), and the stochastic clock is not different from the regular one. 

The situation changes dramatically if the mean of sojourns does not exist, or better to say, it exceeds largely a 
typical time required to diffuse across the physical medium of a finite size. Then, p(n, t) is not a sharp function 
around the mean number (n(t)). This is the case, for example, when possesses a long tail, ip(r) oc (r/r sc ) _1_a 
for t/t sc — > oo with < a < 1. This implies a divergent mean (r) — > oo, and diverging higher moments as well. 
Nevertheless, (n(t)) oc (t/T sc ) a exists for any finite t and it scales sublinearly with the physical time t (see in Appendix 
A). In terms of the one-sided Levy distribution density C a (z), p(n,t) = (t/T sc )£ a (n~ 1 / a t/T sc )/(an 1 / a+1 ). 

Consider now an ensemble of particles. Until time t, each particle has accomplished an individual number of 
intrinsic time periods corresponding to the intrinsic time r{t) which becomes a random variable broadly distributed: 
all the particles have their own history, maintaining individuality and avoiding the fate of self-averaging even in 
the strict limit t — > oo. Only an additional ensemble averaging smears out this principal randomnes o 10 ' 23 i 24 . The 
unbiased diffusion becomes anomalously slow and nonergodic with the spatial variance of a cloud of particles growing 
sublinearly, (5x 2 (t)) cx (n(t)) cx t a . Such a nonergodic approach to subdiffusion seems appropriate for disordered 
solids, e.g. thin amorphous films 2 "! 4 -, with a more recent example provided by Ti02 nanocrystalline electrodes in the 
Gratzel's photovoltaic cell elements^. For charged carriers within such media one can create a potential energy profile 
U(x) by applying an electrical field of spatially distributed fixed charges and a static external electrical field. Then in 
the continuum approximation subdiffusion can be described by the fractional Fokker-Planck equation (FFPE) & 26 i 27 
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where f(x) = —dll(x)/dx is the force, r\ a is the fractional friction coefficient related to the fractional subdiffusion 
coefficient n a by the generalized Einstein relation, r\ a = ksT / ' n a , at temperature T and 

-7 , „ 1 d /"* P(x,t') 
t0 Dl P(x,t) = — -— dt' , ' , (2 

r(i - 7 ) dt J t0 (t - t'p y ' 

is the Riemann-Liouville operator of the fractional derivative 2 ^, where < 7 < 1 and T(x) is the gamma- function. 
The FFPE (fl]) can be derived within the above continuous time random walk (CTRW) framework. It can also be 
written in the form using the Caputo fractional derivative 2 ^ 

«„ W p (m);= ^/'^~^ (3) 

acting on the left hand side, yielding 29 

oU?P( aS ,*) = «a^^W^^Wp( a : 1 t)j =-^, (4) 
in the transport form. Here, /3 — l/(ksT) is inverse temperature and 

J(x,t) = -^e-wl^Wp^i) (5) 

is the subdiffusive flux. It should be emphasized that a non-Markovian Fokker-Planck equation never defines the 
corresponding non-Markovian process completelySi^ 2 -. It allows to find merely the single-time, conditional, and 
double-time probability densities, but never the multi-time probability densities. However, the FFPE dynamics can 
be nicely simulated from the underlying CTRW with the nearest neighbors jumps onl y 29 i 30 . 
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This does not exclude infinite range memory leading to a weak ergodicity breaking-?. 
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A quite different approach to subdiffusion is associated with the fractionai Brownian motion 3 - 3 -. Here, the principai 
issue is the long-range anticorrelations in the particle displacements, positive increments follow with a greater proba- 
bility by negative increments and vice wrs a 34 ' 35 , which can reflect e.g. the phenomenon of viscoelasticity in complex 
glass-forming liquids above, but close to the glass transition. The Brownian particle is temporally trapped in a trap 
(cage effect) by an elastic force with spring constant G(t) which decays to zero in time releasing the particle. Let's 
assume that the motion starts at to, v(t) — x(t) — for t < to- In the linear approximation, on the particle acts 
a viscoelastic force F v _ e \(t) = — J t G(t — t')x(t')dt' , where x(t) is the particle's instant velocity. The first theory 
of viscoelasticity has been proposed by J. Clerk Maxwell 3 ^ in 1867. It corresponds to G{t) exponentially decaying 
in time, G(t) — Go exp(— f/ro), with a relaxation time constant tq. Departing from the phenomenon of elasticity in 
solids Maxwell derived the phenomenon of viscosity in liquids in the limit where the decay of elastic modulus is very 
fast on the time scale of v(t) change, which corresponds to G(t) — 2?7o<5(f), with ry = GqTq being the viscous friction 
coefficient. In the theory of generalized Brownian motion, G(t) is interpreted as the frictional memory kernel 77(f), 
rather than a decaying elastic force constant. Here, one departs from the phenomenon of viscosity and viscous Stokes 
memoryless friction and considers the emerging elasticity in complex fluids or viscoelastic bodies. Both view points 
are essentially equivalent for a positive 77(f) > departing just from different standing points 2 . The anticorrelations 
in the particle's displacements are due to the elastic restoring force component. 

In complex media, the memory function G(t) is better described by a sum of exponentials reflecting a viscoelastic 
response with multiple time scales. Moreover, in 1936 A. Gemant 37 found that some viscoelastic bodies are better 
described by a G{t) relaxing in accordance with a power law, G(t) oc f~ Q , rather than a single-exponential and 
introduced a fractional integro-differential in the viscoelasticity theory. Using the notion of fractional Caputo derivative 
such a visco-elastic force can be short-handed, written as 



F v _ d (t) = -VatoDZxQb) . (6) 

Indeed, such and similar viscoelastic responses are measure d 16 ' 20 ' 38 using the microrheology methods 39 . The Brownian 
motion never stops and the frictional loss of energy must be compensated on average by the energy gain provided 
by a zero-mean random force of environment so that at the thermal equilibrium the equipartition theorem holds, 
in accordance with the classical fluctuation-dissipation theorem. Within the considered model of a linear memory- 
friction such a force must be Gaussian^ - (but not necessarily so beyond the linear friction model). As a result, the 
Brownian motion of a particle of mass m is described by the Fractional Langevin Equation (FLE)i±rii 



mx + Va oD* x{t) = f(x) + £(f) , (7) 
(from now on we fix to = 0) which is a particular case of the celebrated Generalized Langevin Equation (GLE)i£~— 

mx+ [ r](t - t')x(t')dt' = f(x) + £(t) , (8) 



with the memory kernel 77(f) = r\ a t a /r(l — a) and the noise autocorrelation function obeying the fluctuation 
dissipation relation 

mm) = k B T V (\t-t'\). (9) 



Such a GLE can be derived also from a Hamiltonian model for a particle bilinearly coupled with coupling constants 
d to a thermal bath of harmonic oscillators with masses m, and frequencies 0Ji, HB,int(Pi,qi,x) = (1/2) YliiPi l m i + 
rriiUjf[qi — Cix/(miLu 2 )] 2 }. The total effect of the bath oscillators, which are initially canonically distributed with 
Hb, int at temperature T and fixed x = a;(0), is characterized by the bath spectral density 

j( W ) = E£;_^_a( w - Wi ). (io) 

2 * — ' rriiUJi 

i 

The memory kernel is 77(f) = (2/tt) f Q °° J(lj) cos(u)t)du in terms of J(w) and the subdiffusive FLE corresponds to 
a sub-Ohmic, or fracton thermal bath with J(ui) — r] a sin(na/2)uj a — . Without frequency cutoffs such a model 



2 r](t) can also be negative, e.g. accounting for the hydrodynamic memory or in the case of supcrdiffusion. Therefore, the memory-friction 
interpretation is, in fact, more general. 
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presents a clear idealization. There always exists a highest frequency of the thermal bath and this leads to a small 
time regularization of the memory kernel, i.e. a short-time cutoff. Physically, this takes into account the medium's 
granularity beyond the continuum approximation. Moreover, in the case of a finite-size medium there always exists 
also a smallest frequency of the medium's oscillators corresponding to the inverse size of the medium. These facts 
become especially clear if one evaluates the spectral density of low-frequency "fracton" oscillators in proteins, see in 
Refsj^i. This leads also to a cutoff at large times in the memory kernel and the dynamics can be subdiffusive on the 
time scale smaller than the corresponding memory cutoff. The latter one can be but prominently large which makes 
the considered idealization relevant. Important is also the result that an overdamped FLE description of subdiffusion 
can be derived from a broad class of phenomenological continuum elastic models^. 

In the inertialess limit with m —> 0, one can conceive the idea that FFPE ([TJ [3J is the fractional Fokker-Planck 
equation corresponding to the FLE (|7J). This idea is but wrong 4 ^. The non-Markovian Fokker-Planck equation 
(NMFPE) which corresponds to the GLE (with arbitrary kernel)^— and to the FLE, in particular—, is a different 
one. Presently, its explicit form is known only for constant or linear forces f(x)51r— . This resulting NMFPE has the 
form of Fokker-Planck equation with time-dependent kinetic coefficients. This time-dependence is not universal and 
it heavily depends on the form of potential. In turn, the Langevin equation which corresponds to the above FFPE is 
known and it has the form of a Langevin equation which is local in the stochastic time r(t) and describes thus a doubly 
stochastic process^. Here lies also the profound mathematical difference between these two approaches to subdiffusion. 
The physical differences are also immense. In particular, the GLE and FLE approaches are asymptotically mostly 
ergodic as they are not based on the concept of fractal stochastic time with divergent mean period and the mean 
residence time in a finite spatial domain remains finite. Before we discuss the striking differences in more detail, let 
us start from some apparent, but misleading similarities. 



II. FREE SUBDIFFUSION AND CONSTANT BIAS 



Free subdiffusion, as well as diffusion biased by a constant force F can readily be solved in both approaches using 
the method of Laplace-transform. First one finds the Laplace-transform of the mean ensemble-averaged displacement 
(Sx(t)}, and of the position variance (5x 2 (t)) = (x 2 (t)) — (x(t)) 2 , starting from a delta-peaked distribution at x = 
and to = 0. Then, one transforms back to the time domain. This gives 6 

(6x(t)) = n a Ft a /r(l + a) (11) 

and 

(Sx 2 {t)) =2Kat a /T(l + a) (12) 

with the generalized mobility \i a — 1/^q related to the subdiffusion coefficient at F — by the generalized Einstein 
relation /i Q = K a /(ksT). Within the FLE approach these results are valid in the strict inertialess limit m — >• 0. 
Furthermore, the Eq. (| 12[) is still valid then for arbitrary F ^ 0. However, within the FFPE approach the Eq. (fT2]) 
is valid only for F — 0, which is the first striking difference, see also below. Furthermore, both results are also valid 
asymptotically, t — > oc, within the FLE for a finite m/0. 

Generally, the GLE results can be obtained for arbitrary memory kernel r)(t). Assuming the particles being initially 
Maxwellian distributed, i.e. thermalized with thermal r.m.s. velocities 



Vt 



= Vk B T/m, (13) 



one can obtain for the Laplace-transformed stationary velocity (fluctuation) autocorrelation function (VACF) K v (t) = 
(Sv{t + T)Sv(t)}, Sv{t) = v(t) - (v(t)), 

K v (s)= k f. (V (14) 

where fj(s) is the Laplace-transform of r/(t). This is a well-known result which was obtained first by Kubo 45 ' 46 in the 
Fourier space. For the FLE with fj(s) — ?7 a s a_1 it yields by the inversion to the time-domain 41 - 

K v {t) = v 2 T E 2 _ a [-{T/T v ) 2 - a ] (15) 

with t v = (m/r) a ) 1 '( 2 ~ a > being the anomalous velocity relaxation time constant. In ([T5|) . E~ 1 (z) is the Mittag-Leffler 
function, E 1 (z) = X)^Lo z "/r( n 7 + 1)^- For < a < 1, K v (t) is initially positive reflecting ballistic persistence due 
to inertial effects and then becomes negative (anti-persistence due to decaying elastic cage force). In the limit m — > 0, 
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the VACF undergoes a jump starting from Vj, at r = and then becoming negative, K v (r) oc — 1/t 2 ~" for r > 0, 
corresponding to purely anti-persistent motion. The position variance is given by the doubly-integrated VACF. Its 
Laplace-transform therefore reads, 

(fe 2 ^)) = 9r 2fcj?r _. .. . (16) 
v v " s 2 [ms + r](s)} v ' 



Moreover, 



(5x(s)) = F (17) 
s z [ms + rj(s)\ 



for arbitrary kernel, which can also be easily shown from the GLE, and therefore 3 

(Sx(t)) F 
(Sx 2 {t)) 2k B T 



(18) 



for the thermally equilibrium initial preparation. For the FLE with a finite m the inversion of Eq. (|16l) to the time 
domain yields^, 

(5x 2 (t)) = 2v 2 T t 2 E 2 ^ 3 [-(t/T v ) 2 - a ], (19) 

where p(z) = Y^=o zn /^( n 1 + P) 1S the generalized Mittag-Leffler function. One recovers Eq. (JT2J) in the limit 
m — y 0. 

However, for the subdiffusive CTRW and FFPE dynamics the behavior of the ensemble-averaged variance is very 
different from Eq. (fT2"|) under a non-zero bias F ^ 0. Then, the Eq. (JT5J) is not valid anymore. This fact is ultimately 
related to the properties of the stochastic clock. The point is that starting from a CTRW picture it is easy to show (see 
Appendix A) that the growing ensemble-averaged variance (Sx 2 (t)) depends in the asymmetric case (the probabilities 
to jump left and right are different) not only on the mean number (n(t)) of the stochastic clock periods passed, but 
also on their variance (Sn 2 (t)). For a = 1 (regular clock), (Sn 2 (t)}=0. However, for < a < 1, (6n 2 (t)) oc t 2a and 
this dramatically changes the character of anomalous CTRW and FFPE diffusion in the presence of bias. It becomes 
asymptotically (Sx 2 (t)) oc F 2 t 2a , while (Sx(t)} oc Ft a . Notice that for 1/2 < a < 1 the subdiffusion at F = 
transforms into superdiffusion for F ^ 0, i.e. a cloud of particles spreads out anomalously fast relative to its center 
of mass. This yields a remarkable scaling for the ensemble-averaged quantities 

,. (Sx 2 (t)) ,. (6n 2 (t)) 2r 2 (a + l) „ 

im , = im V ; \/J = — — i f - 1 . 20) 

t^oo(Sx(t)) 2 t^oo ( n (t)) 2 r(2a + l) V ' 

This scaling, which was observed first in Refsi^ for a CTRW subdiffusion in the absence of any additional potential 
U(x), has been shown to be universal within the FFPE description also for arbitrary tilted washboard potentials 
and temperatur o 29 ! 30 . Recently, this astounding fact has been related to the universal fluctuations of anomalous 
mobility and weak ergodicity breaking^ 7 . Ultimately, this is just the property of the stochastic clock and it reflects 
the scaling between the variance and the mean number of stochastic periods passed within the external observed time 
t. Surprisingly, the viscoelastic GLE subdiffusion also exhibits a universal asymptotical scaling in tilted washboard 
potentials. In the t — > oo limit it is the same as in Eq. (fl8|) . Astonishingly, it works both for a vanishingly small F, and 
for an arbitrary strong bias. Moreover, both the diffusion and drift in the tilted washboard potentials do not depend 
asymptotically on the amplitude and the form of the periodic potential in the case of GLE subdiffusion and are given 
by Eqs. (fT2j) and (fTTj) . correspondingly. This again is very much different from the FFPE case, where Eq. (p~8|) can 
be used only to calculate the anomalous flux response at a vanishingly small F from the equilibrium (8x 2 (t))p —q at 
F = 0. Also, given (8x(t)) at F one can calculate (Sx 2 (t))F = o using Eq. (fT5)l and the corresponding subdiffusion 
coefficient in periodic potentials in the limit F — > 0, for details see in the work^ 2 - and below. 



For nonequilibrium initial preparations this result holds asymptotically in any asymptotic ergodic case, including the FLE dynamics. 
The relaxation to the asymptotic regime, or aging, can be but very slow 55,56 which is a general feature of subdiffusive GLE dynamics 
also in periodic potentials. 
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III. OTHER SIMILARITIES 

One more similarity emerges for the relaxation of mean fluctuation from equilibrium in harmonic potentials, U (x) = 
kx 2 /2. Then, both the FFPE approach and the FLE approach (in the limit m — > 0) yield the same relaxation 
la w 19 i 26 i 43 , (8x(t)) — (5x(0))E a [— (t/r r ) a ] with the ultraslow position relaxation time constant r r = (rja/k) 1 / 01 . 
Asymptotically, this relaxation follows a power-law, (Sx(t)) oc t~ a . 

The asymptotic distributions of the residence times within a half-infinite spatial domain (or the first return times to 
the origin in the infinite domain) in the case of free subdiffusion are also similar, following the same scaling la w 14 ' 50 i 58 
^(t) cx 1/t 2 ~"/ 2 . However, here the similarities end. The asymptotics for a finite-size domain cannot be same. 
In particular, the mean residence time in any finite-size domain within the subdiffusive GLE description is finite^, 
whereas within the FFPE description is not, except for the case of injection of diffusing particles on the normal 
radiative boundary, where they can be immediately absorbed 14 . Moreover, the GLE (for arbitrary 77(f), including 
FLE) describe a Gaussian process for constant and linear forces f(x) 4 , whereas the FFPE does not correspond to a 
Gaussian process in these cases, see in Ref£. 

IV. DIFFUSION AND TRANSPORT IN WASHBOARD POTENTIALS 

Let us proceed with the case of washboard potentials, where the differences between the two discussed approaches to 
subdiffusion become particularly transparent. We consider the tilted potential U(x) = V(x) — xF, where V(x + L) = 
V(x) is a periodic potential with the spatial period L. 

A. FFPE dynamics 

In this case, one can find exact analytical results for the ensemble-averaged nonlinear mobility (J> a (F) using Eq. (|II[) 
asymptotically also in washboard potentials. First, one finds the exact analytical expression for the ensemble-averaged 
subvelocity v a (F) = fi a (F)F. The FFPE in the form @ is more convenient for this purpose. Indeed, it has the 
form of a fractional-time continuity equation with the flux J(x). For the sake of generality we consider its further 
generalization with a spatially-dependent subdiffusion coefficient K a (x), 

J(x,t) = -K a (x)e-^ u ^ 4- ^ U(x) P(x,t) (21) 
ox 

which is assumed to be periodic with the same period n a (x + L) = n a {x), and the generalized Einstein relation is 
fullfield locally at any x, K a (x) — ksT /r) a (x). We proceed similarly to the case of normal diffusion^— , a = 1. 
A spatial period averaged density P(x,t) = J2k=-M P( x + kL,t)/(2M + 1) should attain a steady-state regime 
(corresponding to a non-equilibrium steady state for F ^ 0) in the limit M — > 00, t — > 00 and that becomes periodic 
with the period L, P s t{x + L) = P s t{x). The corresponding subdiffusive flux J(x), defined with P s t(x), becomes a 
constant J a in the steady state: 

J a = - Ka ( x )e-^^e^P st (x) . (22) 
ax 

Then, the dynamics of the averaged mean displacement follows as 

D?(x(t))=LJ a , (23) 

which can be shown akin to the normal diffusion case^. The appearance of the fractional Caputo time derivative in 
the lhs of Eq. (|23|) is the only mathematical difference as compared with the normal diffusion case. The solution of 
yields for the mean excursion 

(x(f))=«i wb )(F)r/r(l + «), (24) 



4 This is just by the linearity of the transformation from the Gaussian noise to the stochastic process x(t) as described by Eqs. S71 
and JHJ. 
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with wi wb ' ) (F) = LJ a being the subvelocity in the washboard potential. 

One finds J a and v^" h \F) by multiplying Eq. (|22l) with e^ u ^/n a (x) and integrating the result within one spatial 
period. Taking into account the spatial periodicity of V(x) and K a (x) this yields: 

rv+ L P PU(x) 

J a / = -e^y^P st (y + L)+ efOMp^y) 

= (l-e-^ L )e^P st (y). (25) 

Next, multiplying (1251) with e~^ u ^ y \ integrating over y within [0, L], and using the normalization P st (y)dy = 1 
one finds the main result 



v. 



W(F) = { - - '— (26) 



Accordingly, the nonlinear anomalous mobility is ^^^(F) = v^ vb \F)/F. This presents a further generalization 
of the result for subvelocity in Refs . 29 i 30 to a spatially-dependent subdiffusion coefficient n a (x). The subdiffusion 
coefficient in the unbiased washboard potential for F — can also be found using the generalized Einstein relation 
K { a b) (F = 0) = k B T^ h) (F = 0). It reads, 



L 2 

Jo e pV{y >dyJ Tjffdx 



and for k q = const this is the result of the work^ 2 -. For constant n a and a number of different potentials V(x), 
temperatures T and biasing forces F, these two general results were beautifully confirmed by numerical simulations 
of the underlying CTR W 29 ' 30 i 62 on a lattice from which the FFPE in the form (J4j) was derived in the work 2 ^. These 
simulations also confirmed the universality of the scaling relation (|20p within the FFPE approach. Surprisingly, it 
remains invariant also in the presence of a driving which is periodic in time, in the biased case F ^ O^ 3 -, featuring thus 
the universality class of subdiffusion governed by a stochastic clock with divergent mean period and characterized 
by the only parameter a. The above v^ h ^ is the ensemble-averaged result. The subvelocities of individual particles 
remain randomly distributed in the limit t — > oo and they follow a universal subvelocity distribution which reflects 
the distribution of random individual time of travelling particles, as it has been clarified in Ref. 57 . Both the weak 
ergodicity breaking and the universal fluctuations of anomalous mobility within the FFPE approach are ultimately 
related to this remarkable property of the stochastic time. 



B. GLE dynamics in periodic potentials 

The GLE subdiffusion distinctly differs in the physical mechanism and this leads to quite different results for 
washboard potential a 34 ! 64 . First of all, it is asymptotically ergodic and self-averaging over a single trajectory yields a 
quite definite non-random result2i. No additional ensemble averaging is required. Moreover, it turns out that both 
the particle anomalous mobility /Lti wb ' ) and the subdiffusion coefficient k 1 ^^ do not depend asymptotically neither on 
the potential V{x), nor on the bias F being universal and the same as for biased GLE subdiffusion in the absence of 
periodic potential, obeying the generalized Einstein relation. The transition to this asymptotic regime is, however, 
very slow and it strongly depends on the amplitude of the periodic potential Vq and the temperature T. Because 
of this slowness of the transient aging, this asymptotic regime will not necessarily be relevant on a finite time scale 
for anomalous transport in finite-size systems. This is especially so if the periodic potential amplitude exceeds 
the thermal energy by many times. However, this remarkable property features the very mechanism of the GLE 
subdiffusion, which is based on the long-range velocity and displacement correlations and not on diverging mean 
residence time within a potential well, in clear contrast to the CTRW subdiffusion with independent increments. 
It outlines a quite different universality class of subdiffusion. This is the long-range anti-persistence which limits 
asymptotically the GLE subdiffusion and transport processes in the washboard potentials. Since the mean residence 
time in a potential well is finite 3 ^, a coarse graining over the potential period, which makes the sojourns in the 
trapping potential wells irrelevant, becomes asymptotically possible. In fact, upon increasing the potential height 
the escape kinetics out of a potential well (being asymptotically stretched-exponential) becomes ever closer to the 
normal exponential kinetics 3 —, where it becomes described by the non-Markovian rate theor y 48 i 65 . This does not 
mean, however, that the diffusion spreading over many spatial periods becomes normal. As a matter of fact, in the 
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unbiased periodic potentials the diffusion cannot become faster than the free subdiffusion and this is a reason why the 
asymptotic limit of free subdiffusion is attained. A signature of this universality has been revealed theoretically for 
quantum transport in sinusoidal potentials for the case of sub-Ohmic thermal bath which classically corresponds to the 
considered case of fractional sub-diffusive friction. Technically this was done by using two different approaches, one 
perturbative^ and one non-perturbative based on a quantum duality transformation between the quantum dissipative 
washboard dynamics coupled to a sub-Ohmic bath and a quantum dissipative tight-binding dynamics coupled to a 
super-Ohmic bath 4 —. In the quantum case, there are also tunneling processes which are accounted for. Our numerical 
results for the classical Brownian dynamics indicate, however, that this feature is purely classical and, moreover, it 
is universal, i.e. is beyond the particular case of sinusoidal potentials^ -. It is not caused by the quantum-mechanical 
effects. 

Our numerical simulation approach is also insightful and it can be considered as an independent theoretical route 
to model anomalous diffusion and transport processes. The idea is to approximate the non-Markovian GLE dynamics 
with a power-law kernel by a finite- dimensional Markovian dynamics of a sufficiently high dimensionality D ' 6 — . 
Here, "sufficient" means the following: having subdiffusion extending over r time-decades one finds a D-dimcnsional 
Markovian dynamics whose projection on the (x, v) plane approximates the GLE dynamics over the required time 
range within the accuracy of stochastic simulations, as it can be checked for the cases where an exact solution of the 
GLE dynamics is available (free or biased subdiffusion, subdiffusion in harmonic potentials). Increasing D one can 
cover larger r of experimental interest and the embedding dimension D turns out to be finite to arrive at the asymptotic 
results valid for the strict power law kernel. Needless to say that the practically observed cases of anomalous diffusion 
hardly extend over more than 6 time decades (typically several only) which underpins the practical value of our 
approach. 

We expand the power law kernel into a sum of exponentials 

N 

r(i-a) f=( 

obeying a fractal scaling with Vi = &t)/6* , where b > 1 is a scaling parameter, uq > is high-frequency (short-time) 
cutoff corresponding to the fastest time scale tq = 1/vq in the hierarchy of the relaxation time constants, Tj = rob , of 
viscoelastic memory kernel. C a (b) is a numerical constant to provide abest fit to n(t) = r/ a t~ a /T(l — a) in the interval 
[to, Tob 1 ^^ 1 ]. In the theory of anomalous relaxation similar expansions are well-known^2S. In the present context, 
the approach corresponds to an approximation of the fractional Gaussian noise by a sum of uncorrelated Ornstein- 
Uhlenbeck (OU) noises, £(i) = YliLiCi{t)> with autocorrelation functions, (d(t)Cj(t')) = fcsTKjJy ■ exp(— z^|t — t'\). 
This idea is also known in the theory of 1/f noise 71 . For t > r^b 1 ^^ 1 the tail of (|28p is exponential and the diffusion 
becomes normal for t 3> rob . However, by increasing N one can enlarge the corresponding time scale and even 
make it practically irrelevant. The subdiffusion can be modelled in this way over r = iVlog 10 6 — 2 time decades and 
the corresponding embedding dimension, D = N + 2, can be rather small. Such fits are known to exhibit logarithmic 
oscillations superimposed on the power law^. However, their amplitude can be made negligibly small if to choose 
b sufficiently small, e.g. for b = 2 they become already barely detectable. Nevertheless, even the decade scaling 
with b = 10 suffices to arrive at excellent (within the statistical errors of Monte Carlo simulations) approximation of 
the FLE dynamics by a finite- dimensional Markovian dynamics over a huge range of time scales. Weak logarithmic 
sensitivity of r to b and linear dependence on iV allows one to improve the quality of Markovian embedding at a 
moderate computational price. The choice of Markovian embedding which corresponds to l|28p is not uniqu o 68 ' 69 . A 
particular one is the following^: 



/v 



mv 



■hi = -k t v - ViUi + y/2iSikik B T£ t i(t) , (29) 

where fcj = C a {b)rj a vf /T{1 — a) > and £j(t) are independent unbiased white Gaussian noise sources, {£i(t)£j(t')) = 
SijS(t — t'). Indeed, integrating out the auxiliary force variables in Eq. (|29[) it follows that the resulting dynamics is 
equivalent to the GLE (jHJ, (JHJ) with the kernel (f28|) . provided that ttj(0) are unbiased random Gaussian variables with 
variances (uf (0)) = kiksT. The latter condition ensures the stationarity of in the GLE (jHJ, as well as validity 
of the FDR @ for all times. Using different non-thermal preparations of Ui(0) one can study the influence of initial 
non-stationarity of the noise £(t) in the GLE on the Brownian dynamics^ 9 -. In this aspect, our approach is even more 
flexible and more general than the standard GLE approach. 

The auxiliary variables Ui can be interpreted as elastic forces, Ui — —ki{x — Xi), exerted by some overdamped 
particles with positions Xi, which are coupled to the central Brownian particle with elastic spring constants ki and 
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FIG. 2: Anomalous diffusion in the potential U(x) = — Vo sm(2-wx / L) — Fx for various Vo and F at T = 0.1 for a — 0.5. 
Notice an excellent agreement (differences practically cannot be detected in this plot) of simulations with the exact FLE result, 
(8x 2 (t)) = 2v?t 2 E2- a ,3[—{t / T v ) 2 ~ a ], in the absence of periodic potential Vo = 0. Scaling: time in t v = {m/ri a ) l ^ 2 ^ a \ l distance 
in L, energy in m(L/T v ) 2 , force in mL/r 2 and temperature in mL 2 /(r^fc_g). 

are subjected to viscous friction with frictional constants iji — ki/vi = C a (b)rj a T^ a /T(l — a) and the thermal 
random forces of environment. This corresponds to motion of N + 1 particles in a potential U(x, {xi}) = U(x,t) + 
(1/2) J2i=i k%{x ~ Xi) 2 . The Brownian particle is massive (inertial effects are generally included), all other "particles" 
are overdamped (massless, m, — > 0). For example, one can imagine that some coordination spheres of the viscoelastic 
environment stick to the Brownian particle and are co-moving. Their influence can be effectively represented by N 
"quasi-particles" . In this insightful physical interpretation, our embedding scheme is equivalent to: 

N 



= f(x,t) - ^ki{x - Xi) , 



r]i±i = h(x - Xi) + y / 2r/ l k B T£i(t) 



(30) 



It worth to notice that in this approach the mass of the Brownian particle and therefore the inertial effects are 
important. In order to perform an overdamped limit m — > 0, one has to include the viscous frictional force —t/qx 
acting directly on the particle and the corresponding random force. Then, in the limit m — > 0, one obtains 



N 



m x 



f{x,t) - y^^kj(x - Xi) + y / 2r/ k B T^ Q (t) , 



r}i±i = h(x - x t ) + y/2rj l k B T£ i (t) 



(31) 



where (,o(t) is a zero-mean Gaussian random force of unit intensity which is not correlated with the set {£i{t)}. 
However, it was noticed^ that the inertial effects are generally important for the subdiffusive GLE dynamics and 
therefore we take them into account. Of course, here emerges one more difference with the alternative description of 
subdiffusion within the FFPE ©, Q. 

A proper fractal scaling of coefficients ki and r\i with i (see above) allows one to model viscoelastic subdiffusion over 
arbitrary time scales of the experimental interest. One can numerically solve these stochastic differential equations 
(|2T)1) e.g. with a standard stochastic Heun method^ (second order Runge-Kutta method) as done in Refs^"- 64 . An 
example of such simulations is given in Figs. [UGIlfor ol = 0.5, vq — 100, b = 10, C a (b) = 1.3, N ~ 12 and k B T = 0.1. 
The following scaling is used: time in the units of t v , distance in the units of L. All the energy units are then 
scaled in AE = m(L/T v ) 2 and the force units in mL/r 2 . Stochastic Heun method is used to integrate Eq. (|2"9"]l with 



5 This is a natural scaling of the velocity autocorrelation function in time. Other scalings are but also possible 34 ] 64 . They are more 
suitable to consider dynamical regimes close to overdamped. 
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FIG. 3: Anomalous transport for F > 0, see Fig[2]for details. 



a time step At = (1 — 5) • 10 -3 until t max = 2 • 10 5 and n = 10 4 trajectories are used for the ensemble averaging. 
Stochastic numerics are compared against the exact results for the free subdiffusion and for the mean displacement 
under a constant biasing force. The agreement is excellent. The considered particular embedding still works as an 
approximation to the FLE dynamics until t = 10 8 . If one needs to describe subdiffusion on an even longer time scale 
one can increase N. If one needs a better precision of approximation one can make b smaller. Initially all the particles 
are localized at the origin, x = 0, with the velocities thermally distributed at the temperature T . For the time span 
t < t v the motion is always ballistically persistent (superdiffusion) . This reflects the inertia of the Brownian particle. 
It assumes the subdiffusive character for t ^> t v , when the VACF is negative. The presence of a periodic potential 
V(x) — — Vb sin(27ra:/L) dramatically changes both subdiffusion, (x 2 (t)) — (x(t)) 2 , as well as subdiffusive transport, 
(x(t)), on intermediate time scale. However, the long time asymptotics of free or biased subdiffusion are gradually 
attained. The initial behavior still within one potential well remains ballistic. One can conclude that both subdiffusion 
and subdiffusive transport are indeed asymptotically insensitive to the presence of periodic potential within the GLE 
approach. This finding is in a striking contrast with the FFPE approach. However, the transient to this asymptotic 
regime can be very slow, depending on the amplitude of the periodic potential and temperature. 

An interesting phenomenon is also accelerated subdiffusion occurring on an intermediate time scale in tilted wash- 
board potentials, as compared with the free subdiffusion. It can be detected in Fig. [2] for a strong yet subcritical 
bias F = 1 < F cr = 2ttVq/L k 1.2566.... This calls to mind the acceleration of normal diffusion in tilted washboard 
potentials— . However, this accelerated subdiffusive phenomenon occurs only on an intermediate time scale because 
asymptotically the GLE subdiffusion is not sensitive to the presence of the potential. One more interesting effect 
occurs for the initially ballistic transport. It first seems paradoxical that in the trapping potential the initial transport 
becomes faster than in the absence of potential and not vice versa, see in Fig. [3] The result can be understood in 
view of the fact that the minimum of the potential under the strong bias F is essentially displaced in the direction 
of biasing force and the particles are initially accelerated by the additional to F force stemming from the periodic 
potential. 



With this Chapter, we reviewed and scrutinized two different approaches, the FFPE approach and the FLE ap- 
proach, to anomalously slow diffusion and transport in nonlinear force fields with a focus on applications in tilted 
periodic potentials. In spite of some similarities in the case of constant or linear forcings it was shown that the 
nonlinear dynamics radically differ, obeying asymptotically two different universality classes. A first one reflects the 
universal fluctuations of intrinsic time clock and is closely tight to a weak ergodicity breaking. In contrast, within the 
GLE and FLE approach the long-range antipersistence of the velocity and position fluctuations renders the asymptotic 
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dynamics ergodic. One approach seems more appropriate for the disordered solids, or glass-forming liquids below the 
glass-forming transition, as characterized by the nonergodic glass phase. Another one seems more appropriate for the 
regime above but close to the glass transition, or for crowded viscoelastic environments like cytosols in biological cells. 
We have left out further pronounced differences between the FFPE and FLE approaches in the case of time-dependent 
fields ^ 63 ! 74 " — . We are confident that our results not only shed light on the origin of profound differences, but also will 
stimulate a further development of both approaches to subdiffusion, and possibly other interrelationships emerging in 
random potentials. 
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Appendix A: Continuous time random walk and random clock 

Consider a lattice with period a and a particle jumping with probabilities q + and q_, q + +q_ = 1, to the neighboring 
sites after a random clock characterized by the residence time distribution (RTD) iP(t) "ticked" on the next jump. 
Within the physical time interval t there will be a variable random number of intrinsic time periods n(t). The 
probability to make m steps forward and n — m steps backward after n periods is given by the binomial distribution, 
P(m,n) — n\/[m\(n — m)\\qTq t l~ m - Using it one can calculate the first two moments, (x k ) — a fc X)m=o(^ TO — 
n) k P(m, n), k = 1, 2, of the particle displacement after n periods: 

(x(t)) = a(q + - q-)n{t) 
(x 2 (t)) = a 2 [n 2 (t)(q+ - q-) 2 + 4n(t)« 7+ « 7 _] . (Al) 

They are still random quantities because of the randomness of n(t). Each particle has an individual number of periods 
completed until t. For the additional ensemble average one obtains 

((»(*))) ens = a (q+ - q-)(n(t)) ens 
({x 2 (t))) ens = a 2 [(n 2 (t)) cns (q+ - q_) 2 + 4(n(t)) cns q + q_] (A2) 

and for the ensemble-averaged variance (([x — ((x)) cns ] 2 )) cns — ((x 2 )) cns — ((x)) 2 ns 

({Sx 2 (t))) ens = a 2 [{Sn 2 (t)) cns (q+ - q_f + 4(n(t)) cns q+q_] (A3) 

where (Sn 2 (t)) cns — (n 2 (t)) cns — (n(t)) 2 ns is the variance of random periods passed. Notice that ((x 2 )) cns — {{%)) 2 ns 
(([x — (x)] 2 )} cns . Clearly, for a regular clock, {6n 2 (t)) cns = and the corresponding contribution to the ensemble- 
averaged position variance is absent. To simplify the notations, we further denote the ensemble averages as (...) rather 
than ((...)) cns . 

The physical time t can be measured by the sum of independent stochastic periods 7% already completed and one 
not yet completed period t* +1 , t = Y^k=\ T k + T n+ij with n = 0, 1, 2, oo. Therefore, the probability distribution 
p(n,t) to have n time periods within t, 'Y^ =0 p{n,t) = 1, is the n + 1-time convolution of the RTDs ip( T ) ( n times) 
and of the survival probability $(r) = f°° ip(r)dT. Its Laplace-transform reads 

p( n>a ) = izM$( fl )]» (A4) 

in terms of the Laplace-transformed ^{t). Let's consider tp(s) w 1 — (sr sc ) a for sr sc — > 0, where r sc is a time unit of 
measurements. The continuous spatial limit is achieved when a — > 0, r sc — ¥ with K a = a 2 /r^ c being a constant. For 
a finite r sc , considering the scaling limit n — > oo, sr sc —¥ with n(sr sc ) a being finite, one obtains 

p(n,s) = r sc (sT sc ) a " 1 exp[-n(sr sc ) Q ], (A5) 

where n is considered as a continuous variable and 



r(t) = n{t)n 



(A6) 
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is the intrinsic random time. Notice that for a = 1 one finds, p(n,t) = S(n — t/r sc ) by inversion to the time domain. 
That means to say that -r(t) = t is not random. For < a < 1, p(n,t) can be expressed via the one-sided Levy 
distribution density C a {t) whose Laplace transform reads C a {s) — exp(— s a ). Then, all the moments {n k {t)) can be 
easily found from (|A5[) to read 

1(1 + to] 

In spite of the fact that the mean time interval (r) does not exist all the moments of the intrinsic time r(t) are finite. 
This might seem paradoxical. However, the intrinsic time scales with the number of stochastic periods passed and if 
the mean period does not exist the moments of n(t) are nevertheless finite for any finite t. This is because a frequent 
occurrence of very long stochastic time periods within some fixed t implies a smaller value of n(t). In particular, 

(n(t)> = (t/r sc r/T(l + a), 

(n 2 (t)) = 2(</r sc ) 2Q /r(l + 2a), (A8) 



and 

(Sn 2 (t)) _ (5T 2 {t)) _ 2F 2 (l + a) 
{n(t)) 2 ~ (r(i)) 2 ~ r(l + 2a) 



- 1 (A9) 



is the most important property of the stochastic clock. It is primarily responsible for the discussed universality class 
of the CTRW-based subdiffusion associated with the universal fluctuations, and the weak ergodicity breaking. 
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